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Abstract 


Concepts for long-range air travel are characterized by airframe 
designs with long, slender, relatively flexible fuselages. One aspect 
often overlooked is ground-induced vibration of these aircraft. This 
paper presents an analytical and experimental study of reducing 
ground-induced aircraft vibration loads by using an actively con- 
trolled landing gear. A facility has been developed to test various 
active landing gear control concepts and their performance. The facil- 
ity uses a Navy A6-Intruder landing gear fitted with an auxiliary 
hydraulic supply electronically controlled by servo valves. An analyti- 
cal model of the gear is presented, including modifications to actuate 
the gear externally, and test data are used to validate the model. The 
control design is described and closed-loop test and analysis compari- 
sons are presented. 


1. Introduction 

Market analysis of future commercial air travel has identified a significant trans-Pacific market for 
supersonic transports. Fuselage and wing design of these transports is governed primarily by aerody- 
namics and associated structural loads. Placement and sizing of the landing gear occur at a later stage in 
the design process, with design constraints such as braking requirements, turning radius for ground han- 
dling, ground clearance during rotation, and weight distribution as the main criteria (ref. 1 ). 

Long, slender, flexible fuselage configurations, especially those with a long overhang from the nose 
gear to the cockpit, are susceptible to ground-induced vibration problems, particularly those produced 
by operating over long-period, low-amplitude elevation disturbances on runways. Although in-flight 
vibrations are also a concern, the work discussed herein will address the mitigation of vibrations trans- 
mitted from the ground to the aircraft fuselage. The mitigation is accomplished by embedding a control 
system directly into the landing gear. 


This paper presents results from an activity at Langley Research Center that investigated three 
aspects of actively controlled landing gear: analytical modeling, control system design, and experimen- 
tal validation. This work attempts to improve the fidelity of analytical models so that they can be used 
for control design and experimental demonstration of various control philosophies. Development of an 
experimental facility that permits realistic concepts that can be transferred to commercial applications is 
also a goal. 

Development of landing gear analysis dates back to the late 1950’s (ref. 2). Work has included 
numerical simulation techniques and experimental measurements to validate the various computer pro- 
grams. A significant volume of the work available in the literature deals with military aircraft requiring 
accurate prediction of taxi loads over repaired, bomb-damaged runways (refs. 3-6). A computer simula- 
tion program named HAVE BOUNCE (ref. 6) was developed to simulate the dynamic response of mili- 
tary aircraft over bomb-damaged runways. To validate the computer code, model validation was 
performed at the Aircraft Ground-Induced Loads Excitation (AGILE) (ref. 7) test facility at Wright- 
Patterson Air Force Base. Recently, attention has focused on ride quality during taxi, takeoff, and land- 
ing (refs. 8 and 9). A simulation program developed by Stirling Dynamics (refs. 8 and 9) is a good 
example of new simulation capabilities. This particular computer simulation program has been vali- 
dated with data from commercial transports in Europe. 



Because the primary design driver in landing gear design is impact loading, landing gears are typi- 
cally tuned passively for impact loading upon landing. Ross and Edson (refs. 10 and 1 1) are among the 
first to consider an actively controlled landing gear to reduce landing loads. Their work is the basis for 
the actively controlled landing gear concept described in this paper. Ross and Edson demonstrated the 
benefits of using an actively controlled landing gear system to reduce impact loads upon landing and 
while traversing bomb-damaged runways. Work by Freymann (ref. 12) demonstrated analytically and 
experimentally the benefits of actively controlled landing gear in reducing landing loads and vibrations 
under various runway profiles. In reference 1 3, Daniels presented analysis and test results for a Navy 
A6-Intruder landing gear system. This paper discusses an extension of the work in reference 13 to incor- 
porate active controls. An A6-Intruder gear was used in the laboratory because it was readily available. 
Necessary modifications to the gear are described along with the facility used in the experimental vali- 
dation phase. 


2. Symbols 


a l 

lower chamber cross-sectional area, ft 2 


effective orifice area, ft 2 

a r 

snubber annulus area, ft 2 


snubber orifice area, ft 2 

\ 

upper chamber cross-sectional area, ft 2 

C c 

servo valve effective discharge coefficient 

Q 

main orifice discharge coefficient 

C dS 

snubber discharge coefficient 

c t 

tire damping, lbf-s/ft 

C 

arbitrary scalar constant 

d l 

diameter of lower chamber, ft 

D 0 

effective diameter of main orifice, ft 

D op 

orifice plate hole diameter, ft 

D P 

piston diameter, ft 

D pin (X s ) 

metering pin diameter, ft 

Dr 

snubber chamber effective diameter, ft 

D s 

snubber orifice diameter, ft 

D u 

diameter of upper chamber, ft 

d(s ) 

disturbance signal, ft/s 2 
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Ex 

E ow 

E, 

f 

G(s) 

8 

E, 

k(s) 

L 

E u 

l 

M l 

M u 

m{s) 

N 

P High 

P L 

P Low 

Pm 

P s 

P u 

Qc 

Qideal’ Qreal 

Qo 

Q s 


derived quantity defined for equation (A15) 
friction force defined in equation (A28), lbf 
tire force, lbf 
friction force, lbf 

<y 

system transfer function, 1/ft -s 
gravitational acceleration, 32.2 ft/s 
tire stiffness, lbf/ft 

*y 

controller transfer function, ft -s 
lift forces, lbf 

upper chamber total length, ft 
axle off-center distance, ft 
piston and tire mass, Ibm 
upper mass, lbm 
sensor noise, ft/s 2 
normal cylinder force, lbf 
reservoir high pressure, lbf/ft 

'y 

lower chamber pressure, lbf/ft 

'y 

reservoir low pressure, lbf/ft 

2 

nitrogen charge pressure, lbf/ft 

'y 

snubber chamber pressure, lbf/ft 

upper chamber pressure, lbf/ft 

volumetric flow due to servo valve actuation, ft /s 

ideal and actual flow rates, ft 3 /s 

flow rate through orifice plate, ft /s 

flow rate through snubber orifice plate, ft 3 /s 
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Its) 

(s) 

m 

% 

Vm 

V T 


x„ 


Wg 


x c 

x(s) 


y(s) 


p 

Y 

p 

p 

a 


reference signal, ft/s 2 
indicates Laplace transform 
inertial reference ground input, ft 
total volume of hydraulic fluid, ft 3 
initial nitrogen volume, ft 3 
upper chamber total volume, ft 3 
fluid velocity, ft/s 

lower mass inertial reference position, ft 
initial nitrogen level, ft 
piston stroke, ft 

upper mass inertial reference position, ft 

command voltage, V 

error command, ft/s 2 

measured output, ft/s 2 

streamline elevation, ft 

half thickness of lower bearing, ft 

polytropic gas constant 

friction coefficient 

density, lbm/ft 3 

specific weight, Ibf/ft 3 


Superscripts: 

i indicates inflow 

o indicates outflow 

Subscripts: 


1,2 points along streamline 

L refers to lower mass 
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u 


refers to upper mass 


(-•) refers to second derivative with respect to time 

(•) refers to first derivative with respect to time 

3. Analytical Model 

To extend the work by Ross and Edson (ref. 10), this research discusses an independent develop- 
ment of a main landing gear mathematical model. The nonlinear equations of motion were developed 
for a telescoping main gear modified with an external hydraulic system for actuation and control of the 
gear. Specific details of the landing gear were taken from technical drawings supplied by the Grumman 
Company. 

Figure 1 shows a schematic of a landing gear used in the development of the equations of motion. 
This schematic is representative of a general telescoping main landing gear. The model includes the 
aerodynamic lift L on the airplane, the mass of the airplane’s fuselage, lumped with the mass of the 
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main cylinder a s M lo and the mass of the piston lumped with the mass of the tire as M L . The inertial 
position of the upper mass is X wg with zero value when the gear is fully extended and the tire is just 
touching the ground. From this same configuration, X a is the position of the lower mass taken as zero at 
the axle of the tire. When the gear is compressed, X a measures the deflection of the tire to an inertial ref- 
erence ground input U(t). Part of the upper cylinder chamber is filled with compressed nitrogen to pro- 
vide the system with a spring. The cross-sectional area of the upper chamber is denoted by A u and the 
corresponding pressure is P u . Likewise, the lower chamber has the cross-sectional area denoted A L and 
a corresponding pressure P Hydraulic fluid moves between the upper and lower chambers through an 
orifice plate with a hole of diameter D op . A tapered pin attached to the piston, known as a metering pin, 
is used to obstruct the flow and effectively vary the orifice diameter as the pin moves through the ori- 
fice. The pin diameter is a function of X s and is denoted as D pin (X s ). Hydraulic fluid reaches the snubber 
chamber through several orifices of diameter d y In the snubber chamber, the annulus area is denoted by 
Ar and the pressure is P s . The diameter of the piston is D p . Some of the actual dimensions are defined in 
table 1 . The figure denotes entry-exit ports in the upper and lower chambers for the exchange of hydrau- 
lic fluid used by the active control system. Tire spring and damping coefficients are denoted by K 
and C t . 


Table 1 . Properties and Dimensions for the Navy A6-Intruder 


Landing gear properties 

Dimensions 

D u 

0.1524 m, (6 in.) 

d l 

0.1524 m, (6 in.) 

Ds 

0.0041 m, (0.1614 in.) 

D P 

0.1397 m, (5.5 in.) 

D 0 

0.02859 m, (1.125 in.) 

L 

0.3832 m, (15.087 in.) 

M l 

145.1 kg, (319.22 lbm) 

M u 

4832.7 kg, (10609.94 lbm) 

Y 

1.1 

P 

5 cm, (2 in.) 

P 

912 kg/m 3 , (1.22 x 10 8 lbm/in 3 ) 


Figure 2 shows the forces acting on the upper mass. Balancing the forces acting on the upper mass 
yields the following equation: 

M u^ng = M u g-L-P u A o -P L (A L ~ A o ) + P s A R -f = F l -f (1) 

where F j is a newly defined term in equation (1), g is the gravitational acceleration, and/is the friction 
force between the piston and the cylinder wall. All other terms were described previously. This equation 
assumes that the hydraulic fluid pressure in the upper cylinder is identical to the nitrogen pressure. Also, 
m this development, the variable A 0 , the main orifice area, shows that the metering pin is included; that 
is, it is a variable cross-sectional area depending on the stroke of the piston. 

Figure 3 shows the forces acting on the piston. By adding the forces on the lower mass (piston), the 
force balance equation is 

M L X a = M L g + P l (A l -A s )-P s (A R -A s )-F t + f = F 2 + f (2) 
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Figure 2. Schematic of upper mass and main cylinder. 



+ Ft 

Figure 3. Schematic of lower mass. 


where F 2 is newly defined in equation (2); F t is the force that is transmitted through the tire from the 
ground and has the form 

F, = K f (X a + U) + C f (X a + U) 

where the tire force is defined as a linear function of tire stiffness and damping. The tire stiffness and 
damping coefficients are obtained by linearizing the behavior of the tire about its nominal operating 
point. Because all pressures are functions of stroke, a more convenient coordinate to use is stroke. 
Defining the stroke coordinate as X s = X wg - X a , equations (1) and (2) can be written as 
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( 3 ) 


The discussion so far relates forces F ] and F 2 to the corresponding pressures. Chamber pressures must 
now be related to displacement and velocity of the gear. Appendix A discusses the relationship of pres- 
sure to piston motion as well as pressure control using external tanks. Gear actuation using external 
tanks enables the design of control strategies to reduce unwanted motion due to ground inputs. In the 
following, a description of a general approach for control design is presented. 


4. Control System Design 


To control the motion of the landing gear, hydraulic fluid from auxiliary tanks is used in conjunc- 
tion with electronically controlled valves to actuate the gear. The goal for control design is to minimize 
disturbance propagation from the ground into the fuselage. To aid the discussion on control design 
methodology, consider a linearized representation of the landing gear and servo valves that are trans- 
formed by using Laplace’s transform into G(s). Using feedback control, as indicated in figure 4, one can 
design a controller £(s) to command the servo valves. 

Define r(s ) as an arbitrary input reference signal, d(s) as an unknown external disturbance, y(.v) as 
the controlled response, and m(s ) as sensor noise. After some block diagram manipulation, the con- 
trolled response is given by 


y(s) = [/ + G(j)&(j)] 1 [</(.? ) + G(s)k(s)r(s) - G(s)A;(s)m(j)] (4) 

The factor I + G(s)k(s) is the return difference and multiplies every term in the right-hand side of 
the equation. To minimize the effects of the disturbance d(s) on the response, the factor multiplying the 
disturbance term d(s) must be made small; that is, the return difference must be large, G (s)k{ s ) » 1 , 
in the frequency range of interest. Because G(s) is fixed, the control designer’s task is to maximize the 
return difference value while maintaining the stability of the system. To ensure a stable design, the 
Nyquist criterion is used for this single-input, single-output problem. Because the landing gear behavior 
is highly nonlinear, one must examine bounds of variations in the system dynamics to ensure a stable 
design. The Nyquist criterion was computed experimentally to assess stability and gain margins of the 
design. These techniques provide tremendous insight into design philosophy and stability analysis, 
although their application to nonlinear systems is limited. 


5. Experimental Facility 


Figure 5 shows a Navy A6-Intruder main landing gear installed underneath a drop carriage in the 
standard vertical position. A connecting plate was fabricated to allow for normal mounting of the gear to 
the plate, and the plate was then rigidly connected to the drop carriage. The drop carriage is a truss 
structure that weighs about 4.5 tons and allows unrestrained vertical motion. The drop carriage rests on 
the landing gear. This mass simulates the rigid portion of the aircraft mass carried by the gear. Once the 
gear is loaded, a shaker table is used to input forces into the gear. Hydraulic lift cylinders, powered by a 
hydraulic pump,- are used to lift the drop carriage and unload the gear. Once the gear has been lifted, the 
ability exists to lock the gear in that position with hydraulic valves. 
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The hydraulic shaker table was built specifically for the task of actuating landing gears. The specifi- 
cations included the capability to perform a 1-in. step bump in 2 ms while bearing 12000 lb mass. Input 
waveforms such as 1 — cos(x), sin(x), trapezoidal with user-selected rise time, and a sawtooth waveform 
are all accurately reproduced by the shaker table. General profiles using runway elevation versus time 
data are also reproduced well for low frequencies. The shaker table is capable of applying dynamic 
forces up to 12 000 lbf on the test mass and allows actuator movement of 6 in. 

The landing gear was modified in a number of ways. Two electrohydraulic servo valves were 
attached to the outside of the landing gear on flat areas that had been machined on the outer cylinder of 
the gear. One valve was located above the orifice plate of the gear (in the upper chamber), and the other 
valve was located below the orifice plate (in the lower chamber). Holes were machined into the gear so 
that the valves could transfer pressurized hydraulic fluid either into or out of the desired chamber. Both 
valves were designed to have flow rates of at least 26 gal/min at 600 psi with a response approaching 
100 Hz. A high-pressure accumulator was mounted on the upper mass (drop carriage) and was kept 
charged to a pressure approximately twice that of the static-loaded charge pressure in the landing gear. 
A low-pressure accumulator was also installed so that when desired, pressurized hydraulic fluid in the 
landing gear could be directed there, reducing the transient back pressure that would tend to restrict the 
outward flow of hydraulic fluid. The low-pressure accumulator was maintained at essentially atmo- 
spheric pressure. Ultimately, the low-pressure accumulator was attached to an atmospheric pressure res- 
ervoir where the pump used to supply the high-pressure accumulator was located. The system was thus 
pressure balanced evenly around the nominal static, loaded charge pressure of the landing gear, permit- 
ting roughly equal flow rates into or out of the gear at similar servo command levels. 

The piston head of the landing gear was also modified. Normally, such a landing gear has a snubber 
chamber that is designed to limit the speed of piston extension to prevent a significant “bottoming out” 
shock on the gear components such as might occur after a catapult during an aircraft carrier launch. 
Thus, normally the hydraulic damping characteristics of the gear vary depending on the direction of pis- 
ton travel. In this experiment, it was desirable to remove the “snubber” effect so that the damping 
behavior was more even in both directions. To that end, a ring mounted directly under the piston head, 
which normally acts as a directional valve and restricts hydraulic fluid motion in one direction, was 
modified by drilling additional holes to provide equal flow past the ring, regardless of the direction of 
hydraulic fluid motion. These changes were accurately reflected in the modeling of the gear for analyti- 
cal purposes. 


The top of the landing gear was modified slightly to accept a high-strength site glass. This site glass 
allowed a visual indication of the proper servicing level of hydraulic fluid, prior to being pressurized 
with nitrogen, and saved significant time in pretest operations. 

The gear was instrumented to provide the necessary information for model validation. There are 
two accelerometers, one placed at the upper mass and the second at the lower mass. Two potentiometers 
were also used, one to locate the upper mass with respect to a fixed position on the carriage and one to 
measure the relative position between the upper and lower masses of the gear. Two pressure transducers 
were used to verify some basic model assumptions, mainly that the hydraulic fluid and the gas do not 
mix to any significant degree after initial shaking. One pressure transducer was located just outside the 
charge port of the upper cylinder, and the other is embedded in the piston head. By using a strain gauge 
on the wheel axle, strut vertical load was inferred by measuring bending moments induced by the tire. 
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6. Test Results and Model Validation 


The following section discusses experimental results from tests conducted on the landing gear sys- 
tem. First, the servo loop dynamics and electronics are characterized and compared with the analytical 
model. Second, the simulation model, which was constructed using a commercially available software, 
is described. Finally, test results for various open-loop and closed-loop cases are presented. 

6.1. Servo Loop Dynamic Characterization 

Figure 6 presents a plot of flow rate as a function of servo command. These data were measured by 
removing nitrogen from the unrestrained landing gear and computing flow rates by measuring piston 
stroke rates as a result of discrete servo commands. Tests were then conducted with the piston restrained 
from moving to characterize the servo loop dynamics with minimum interference from piston motion. 

The slope of the measured flow rate versus command gives the product C c JKp in equations (A 16) 

and (A 17) where A P is the pressure difference between the supply or return and the strut internal pres- 
sure. By using these results, the servo effective discharge coefficient was calculated to be 

C c = 1.0765 x 10~ 6 . To compare simulated chamber pressures to test, a test was conducted using a 

sinusoidal sweep from 0.5 Hz to 10 Hz. Input voltages from the test were input into the simulation, and 
the computed frequency response for upper chamber pressure to servo command is shown in figure 7. 
Test results are depicted with a solid line and simulation with a dashed line. Lower chamber tests results 
(not shown) are similar but with slightly more phase delay between commands and internal pressure 
variations. The initial pressure in the chamber was recorded as 350 psi, the initial stroke was 10.3 in., 
and the high-pressure accumulator pressure was 750 psi. The nitrogen level was estimated to be 4.78 in. 
Simulation results were generated by using Simulink/Matlab (ref. 14) computer simulation code (The 
Mathworks, Inc., Natick, MA). The upper curve in figure 7 shows the magnitude ratio of upper pressure 
to input voltage as a function of frequency, whereas the lower curve shows a phase comparison. This 
transfer function represents the servo valve hydraulic system response at the conditions mentioned 
previously. 



Figure 6. Flow rate as function of servo commands with unrestrained landing gear and no nitrogen. 
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Figure 7. Upper chamber pressure to servo command transfer function. 

6.2. Landing Gear Dynamic Characterization and Model Validation 

Runway elevations and servo command voltages were the two inputs used to characterize the 
landing gear. Because the system is highly nonlinear, sine sweeps were the main form of excitation. 
Time simulations were performed with Simulink. Because the model is nonlinear, initial conditions for 
the different parameters have to be set properly or time integration will fail. Conditions such as upper 
mass position and velocity, piston stroke, stroke rate, upper chamber pressure, and nitrogen level must 
all be specified. In the initial design, two sensors were used to control the motion of the gear, piston 
stroke, and upper mass acceleration. Because the axle load signal from the strain gauge is proportional 
to the upper mass acceleration, the strain gauge output was used for later tests. The axle-mounted strain 
gauges had the additional benefit of being relatively “quiet” and avoided the more dynamic nature of 
acceleration measurements on stiff structures such as those observed with the upper mass 
accelerometer. 

The controller used for all the closed-loop tests was synthesized with the aid of an experimentally 
determined Nyquist diagram. To compensate for phase lag of the servo valves and hydraulic system, a 
lead-lag compensator was used to add about 10° of lead at 1 .5 Hz. Direct axle load feedback with a loop 
gain of 1 V/6731 lbf was used for all the closed-loop test results shown. 

Shaker head position, servo-input command, piston position, upper mass position, internal pres- 
sures, and acceleration responses are compared to simulation results in figures 8 and 9. Solid lines cor- 
respond to test results and dashed lines are simulation results. The input runway elevation is a sinusoid 
with amplitude of 0.75 in. at a frequency of 1 .5 Hz. Piston position feedback is always used to maintain 
interchamber leakage through the servo valve to prevent depleting the hydraulic fluid in the strut. 
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Figure 8. Controlled test and analysis results due to periodic disturbance at 1.5 Hz. 
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Axle load feedback is used as a measure of the upper mass motion. This control loop is toggled on and 
off during an experiment. Data shown in figures 8 and 9 had the axle load feedback loop turned on after 
14 s. The upper mass position is reduced to 25 percent of the uncontrolled position after the axle load 
feedback loop is turned on. Drift after the initiation of control in the stroke and upper mass position his- 
tories in figure 9 could be attributed, in part, to a continuous decrease in the control system hydraulic 
supply pressure. All simulation results assume a constant control system hydraulic supply pressure. To 
experimentally minimize the effect of reduced control system hydraulic supply pressure, long duration 
tests were interrupted periodically to allow for the recovery of system hydraulic pressure. Discrepancies 
in stroke levels between test and simulation are not well understood. 


Friction played a key role in unrestrained tests performed with this test bed. To illustrate the prob- 
lem, figure 10 shows a frequency response function of the upper accelerometer-to-servo command. 
Note that the landing gear locks above 0.7 Hz due to friction. Control authority is lost beyond 0.7 Hz 
due to high friction levels, about 2000 lbf statically and 400 lbf dynamically. This static friction level 
causes a condition in which pressure versus stroke equilibrium can be in error by as much as ±45 psi. 
Also important is the use of nitrogen in the upper chamber. Nitrogen serves as a soft cushion for load 
transfer through the gear. In the absence of nitrogen, the strut is full of hydraulic fluid, which is incom- 
pressible, and therefore small amounts of hydraulic fluid into or out of the strut cause large changes in 
the internal pressures. Because the servo hydraulic system has a limited supply of external hydraulic 
fluid, the absence of nitrogen allows for longer test time and higher forces in the system but requires 
higher pressures for the external supply. 

Controlled tests like the one shown in figures 8 and 9 can be performed only at discrete frequencies 
with the capabilities of the existing hydraulic system. To test the frequency range between 0.1 Hz to 
4 Hz, the spectrum analyzer was set for a sine sweep, and the test was conducted over a long period of 
time, with periodic stops to allow for the hydraulic system to be resupplied. Open and closed-loop 
results from this test are shown in figure 1 1 . Note that in this test the shaker table was used as the input 
disturbance and provided enough energy to prevent the system from locking up below 1.2 Hz. Feedback 
from position and axle load signals were used in the control system to attenuate responses between 1 .5 
and 3.5 Hz. The maximum amplitude reduction is a factor of 4.4 at 1 .4 Hz, with reductions beyond 
3.5 Hz of about 20 percent. Using the strain gauge sensor to measure axle load provides a cleaner signal 
for feedback and reduces the risk of high-frequency instabilities in the feedback loop. Note that in the 
ideal case with the control system fully charged, with a gain optimized for a single frequency, and by 
controlling the system at its natural frequency, amplitude reductions of a factor of 10 have been 
observed. 
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Figure 1 0. Frequency response from upper servo command to upper accelerometer. 



•5 1.0 1.5 2.0 2.5 3.0 3.5 4 . 

Frequency, Hz 

Figure 1 1 . Upper acceleration to shaker position transfer function. 






7. Summary 


Equations of motion for a telescoping landing gear system have been developed, incorporating an 
external servohydraulic system that allows for landing gear actuation. The electronic servohydraulic 
system model combined the electronics and hydraulic dynamics in one relatively simple formulation. 
This formulation uses an exact internal pressure expression requiring the solution of a multivalue alge- 
braic equation at each time step. At equilibrium, pressures within the strut are balanced at a given stroke 
position and command flow into or out of the chamber. 

A number of aspects of actively controlled landing gear design have been demonstrated in this 
study. Fuselage vibration reduction levels, by a factor of 4, have been demonstrated, along with some of 
the fundamental limitations of implementing such systems in landing gear design. High friction levels 
hindered our ability to achieve higher performance without a major redesign of the gear. However, even 
modest vibration reductions may translate into reductions in landing gear loads and therefore aircraft 
structural weight. 
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Appendix A 

Development of Equations of Motion 

Fundamental relationships between forces acting on the landing gear and the fluid pressures inside 
the landing gear are developed in the following sections. 


Pneumatic Pressure Equation 


The pressure terms that contributed to the forces F, and F 2 in equation (3) need to be related to 
position X w? and X s or their derivatives. Nitrogen pressure changes have been described (ref. 13) by 
the polytropic gas law for a closed system as 


P=P 

u m 






(Al) 


where V° ni is the initial nitrogen volume, V ni is the actual nitrogen volume, P ni is the initial nitrogen 
charge pressure, and y is the polytropic gas constant. This representation of the pressure change is 
assumed to happen as a quasi-equilibrium process. In most situations the polytropic gas constant is actu- 
ally not constant and is usually calculated from pressure-stroke data. However, reference 13 showed the 
experimental estimation of a single value for the polytropic constant. Equation (Al) is defined in such a 
manner that P will become very large for small values of V nj ; that is, the gear is nearly collapsed. 
This equation is a suitable representation of the process, with only the polytropic gas constant y as an 
unknown. 


Before developing expressions for hydraulic fluid motion between chambers, a fundamental 
assumption is that the upper chamber hydraulic fluid pressure equals the nitrogen pressure. An expres- 
sion for the nitrogen volume V nj is now sought. The total volume available in the system is defined by 
using the dimensions shown in figure Al as 


V T = A L( L u- X J + A r X s (A2) 

assuming A L = A u . This volume houses both hydraulic fluid and nitrogen in the system. Since the 
nitrogen level is given by X n -, with volume computed as V n - = A j4 X n - , the corresponding volume of 
hydraulic fluid is given by 


V U , = V T- V ni^ A L L u-( A L-\'»‘,-A L X n 


(A3) 


and the rate of change of the volume of hydraulic fluid is 


dV,. 

-Jr = < A L- A r) X s~ A L X » 


(A4) 


In a closed system, the amount of hydraulic fluid V n is constant; therefore, the rate of change is 
zero. Using this fact with equation (A4) yields an expression for the nitrogen level as 


X ni = 


( A L~ A r> 


(A5) 
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Figure A 1 . Hydraulic fluid-nitrogen volumetric diagram. 


If the system is not closed; for example, if hydraulic fluid is entering or leaving the chamber at a rate 
Q , the rate of change of the volume of hydraulic fluid is given as 


^ = Q c = ~(A L -A R )X s -A L X ni 


dt 

Solving this equation for the nitrogen level yields 

v _ Wt-V 

X n i 


Qr 


Xs ~ A~ 

a l 


(A6) 


(A7) 


This differential equation, along with the initial value for the nitrogen level X ° m : , describes changes in 
nitrogen level as a function of time. Canceling the area terms in equation (Al), the upper chamber pres- 
sure is now written as 




P = P . 

u ni 


( Xu ^ 

V X niJ 


(A8) 


Note that equations (A7) and (A8) relate the stroke to volumetric changes in hydraulic fluid. We will 
now develop the pressure equations for the system. 


Lower Chamber Hydraulic Pressure Equation 

Hydraulic fluid pressure in the lower chamber and the snubber chamber is related to the flow rates 
of the hydraulic fluid into and out of those regions. The volumetric flow rates through the orifice plate 
hole O and the snubber orifices Q can be determined by combining the continuity equation and 
Bernoulli’s equation for fluids. Flow is always from high pressure to low pressure. Bernoulli’s equation 

p 1 2 

for an incompressible fluid states that along a streamline - + — v +z = c, where P is the pressure at 

some point, g is the gravitational acceleration, v is the velocity of the flow, d is the specific weight of 
the fluid, z is the height difference from some zero reference, and c is a constant. This equation assumes 
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that the viscous effects within the fluid are negligible, the flow is steady and incompressible, and only 
points along a streamline were considered. Examine points 1 and 2 along a streamline; the mass conti- 
nuity for incompressible fluid states that the flow rates A,v, = A 2 v 2 . Assuming that the pressure at 

point 2 is P 2 > P , and that the flow is from point 2 to 1 , one can solve for the velocity v , from the con- 
tinuity equation as 


where the area in the continuity equation is now expressed in terms of diameters. Substituting this 
velocity into the Bernoulli equation yields 


v 


1 


2(^2- /> l) 

P 

l'-21 


D 2) 


(A9) 


The ideal volumetric flow rate (Q i dea i) for an incompressible fluid can be expressed as Q u j ea i = A r. In a 
realistic flow situation though, there is a flow rate loss due to flow restrictions in the orifice. This loss is 
empirically quantified by a discharge coefficient C d , which represents the percentage of the ideal flow 
that actually occurs. This coefficient, when multiplied by the ideal flow, yields 

Q real = C d Q ideal = AC d v (MO) 

Substituting the velocity in equation (A9) into equation (A 10) 


Q 


real 



J p 2- p vf° rp 2 >p i 


(All) 


For the landing gear shown in figure 1, there are three flows that are of concern: the flow through 
the orifice plate Q a , flow in and out of the snubber chamber Q s , and if the gear is actuated externally, 
there is an additional flow Q c . The Q c is supplied by an external high- and low-pressure reservoir, which 
adds or removes hydraulic fluid from the system. 

Defining a control volume as shown by the dashed line in figure A2, compression occurs when 
> 0 , and extension occurs when X s < 0 . Flow is assumed to be positive leaving the control volume 
and negative entering it. For an incompressible fluid, continuity yields 

Q c + Q 0 + Q s = a l x s (A 12) 

Equation (A12) defines the general form of the continuity equation. Because flow directions have to be 
taken into account explicitly to determine pressures in the various chambers, one needs to consider all 
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chamber 


— 

X, 
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Figure A2. Control volume between piston and orifice plate. 


possibilities. Flow in or out of the snubber chamber always equals the volume displaced during expan- 
sion or compression Q s = A R X S . Substituting this value into equation (A 12) yields 


Q c + Qo = ( a l -A r )X s (A13) 

Substituting the appropriate pressures, areas, and diameters into equation (All), the flow rate through 
the orifice plate during the compression X s > 0 can be written as 


Q = A C 

^ O n i 


O^d 



(d y 

4 

1 - 

o 



D r 



\ L ) 



T'V'V r L > p u 


( A ! 4) 


where D q is the effective diameter of the main orifice, D L is the diameter of the lower chamber, C d is 
the discharge coefficient of the main orifice, and A q is the effective area of the main orifice. To sim- 
plify equation (A14), the nonpressure terms were grouped as 


E, = AC, 

1 o a 


(D\* 


1 - 


\°L) 


The flow through the orifice is now written as 


Qo = £ l V'W. P L> P u 


(A15) 
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Obviously, if the upper chamber pressure is higher than the lower chamber pressure, the order of the 
pressure terms must be reversed. 

To actuate the piston, one can add or remove hydraulic fluid from the landing gear. Considering a 
system with infinitely large high- and low-pressure reservoirs, depicted in figure A2, an external flow 
into the lower chamber can be written as 


Ac = - C c*J > w - P L *c<° < A16 > 

and flow out of the lower chamber is given by 

Q° = CScFFriZ, *c>° < A17 > 

where x c is a command variable, and the parameter C c combines empirically determined values for the 
discharge coefficient, orifice area, and flow gain per unit command input. Superscripts i or o are used to 
distinguish inflow from outflow. To evaluate the pressures, one needs to examine the different flow 
directions in the control volume. 

Figure A3 shows three potential flow directions for Q c and Q 0 for the compression case. The cor- 
responding continuity equations for X s > 0 are 

Q 0 “ Qc = J^L ” P u + C c x cJP High ~ ^ L ~ ( A l ~ x c < Q 0 > 0 
e„ + e” = e , + c c x c jr L ^Z = (a l -a r )x s -, * c >o,a 0 >o (aw 
H2„ + G° = -E, Jp^T L + C'X'Jr’L-Pur = Mj. - a r )X,-, * e >0,Q o <0 

Similarly, the extension case has three potential flow directions shown in figure A4 with corresponding 
flow equations for X s < 0 : 

-Qo-Q> E 1 V'Wl * C c x cJ^H^h~^L = < A L - A *)*,-■ *c<°-Qo< 0 

-e„ + q° = - A » A; < A,9 > 

e„ + Q[ = + C c x 'Jr Htgk -p L = - a k ,x„ x c <0,q o >0 
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Figure A3. Potential flow directions for compression case X s > 0 . 
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Figure A4. Potential flow directions for extension case X s < 0 . 


Equations (A8), (A 18), and (A 19) need to be solved simultaneously to determine the pressures P L and 
P u . Note that equations (A7) and (A8) couple the upper pressure to the flow rate Q c . Appendix B 
shows an explicit solution for the lower pressure P L , obtained by using the Mathematica Symbolic 
manipulation program (The Mathematica Book, third ed., Wolfram Media, Cambridge University 
Press), reference 15. 

A case not considered thus far is when the system is locked due to friction; that is, X s = 0 (lower 
and upper chambers move together). When no external hydraulic fluid source is in use, x c = 0 ; there 
is no hydraulic fluid transfer between chambers. However, when an external source is in use, hydraulic 
fluid transfer is determined by continuity Q 0 = Q . Two hydraulic fluid flow cases are possible 
depending upon the control variable x c value; they are 


-<2o = " £ 1 V'W, = Q'c = • *c < »• Qo < 0 

(A20) 

Qc - Q° - CS'JrL-rL., ■ *.><>. e„>o 


Solving equation (A20) for the lower chamber pressure yields 


/ E ^ 


P L = 


C c x c) 


P u + P High 


( £, y 


for x < 0 

c 


+ 1 


P L = 


f £, 


P +P, 

u Low 


( £, r 

\~ x cj 


for x c > 0 


+ 1 


(A21) 


Most of the pressures needed to solve equation (3) have been developed, except for the snubber cham- 
ber pressure. In the following section, these expressions are developed. 


23 



Snubber Chamber Hydraulic Pressure Equation 

The analysis in this section is similar to the lower chamber analysis just described. To develop pres- 
sure relationships for the snubber chamber, consider a control volume as shown by the dashed line in 
figure A5. The variables A R and D R in figure A5 are the snubber chamber annular area and effective 
diameter, respectively; P is the pressure in the snubber chamber, and D s is the diameter of the snub- 
ber orifices. 

In general, the continuity equation relates the flow in or out of the snubber to the volumetric change 
as follows: 


Qs = A R*s 


(A22) 


Flow into the control volume of the snubber area during a compression mode X 5 > 0 is given by 


Qs = A'C 



(D\ A 


p l> p s 


K D Lj 


= p lJ P ^ rP s 


(A23) 


where C ds is the discharge coefficient of the snubber orifices, and A s is the effective area of the snub- 
ber orifice. For the extension mode X s < 0 , the flow is into the control volume and equation (A23) 
yields 


Q, = 


~ A s C dS 



(D 

4" 

1 - 

s 



d r, 

\ K J 



P . >P L 


- - E iJ P s - P L 


(A24) 



Figure A5. Control volume for snubber chamber. 
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where D R is the effective diameter of the annular snubber chamber. Often the snubber orifice diameter 
D s is different in extension and compression. For those cases, the appropriate diameters need to be 
used in this equation for each case. Equations (A23) and (A24) can be used to determine the snubber 
pressure, given a flow rate Q s . Substituting the flow rate Q s in equations (A23) and (A24) into equa- 
tion (A22) yields 


E 2> 

lr L - F , - v. 

for 

X,>0 

~ E 3, 

A R*. 

for 

X,<0 


(A25) 


Solving for the snubber pressures 


P s = P L~ 


(A*f 


x] for X s > 0 


P=P,+ 


(Ar) 2 -.. 2 


l E, 

V 3 j 


X s for X s < 0 


(A26) 


Equations (A8), (A18), (A 19), and (A26) provide the fundamental relationships between the pressures, 
displacements, and commanded input x . To evaluate the time responses for all these quantities, equa- 
tions (3) and (A7) must be integrated in time with appropriate initial conditions for all variables. Fric- 
tion models will now be described. 


Friction Models 

The only unknown term left in equation (3) is friction. Friction in this landing gear comes mainly 
from two sources: friction due to tightness of the seal and friction due to the offset wheel (moment). The 
seal friction has a maximum value when the system is locked and decreases as the system begins to 
move. The functional relationship between frictional force level and velocity is determined through test- 
ing. The friction due to the offset wheel is the result of the moment produced by the nonaxially loaded 
piston within the cylinder. 

Figure A6 shows that the force between the piston head and the cylinder A is a result of the tire 
force F t applied at a distance / from the centerline of the piston. The frictional force due to the offset 
wheel is 


F ow = \ IN (A27) 

where N is the normal force on the cylinder wall resisting motion of the piston head, and (-1 is the coeffi- 
cient of friction between the two parts. The normal force N is computed by adding the moments about 
point O; P is defined as one-half the thickness of the lower bearing. Solving for the normal force N and 
substituting into equation (A27) yields 


F = n 

ow 


r IF f ^ 

vVP; 


(A28) 
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Cylinder 



Figure A6. Schematic of landing gear friction model development. 


The total friction force in the landing gear/ is given by 

f = F seal +F o W (A29) 

Friction effects were included in the numerical simulations by using the approach described in 
reference 16. 

Note that this model takes into account only vertical loads on the strut. Furthermore, the tire is mod- 
eled as a nonlinear spring and damper combination and does not take into account radial stiffening due 
to centripetal forces. Also, all structural members were assumed to be rigid, each having only a vertical 
degree of freedom. These assumptions are adequate for taxiing over runway profiles and for landing 
impact (spin-up drag on the tire does not significantly affect the vertical loads on the strut). Any braking 
or turning maneuvers were not covered in the development. The equations developed here are the basis 
for a “rollout” simulation. 
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Appendix B 

Analytical Solution of Pressure Equation 

All possible sign combinations in equations (A 18) and (A 19) can be considered by using a generic 
form such as 


~ P u + C 2 J P L~~ P r “ C 3 

(Bl) 

- p L + c 2 J P L ~ P r = c 3 

(B2) 

- p L + C 2 J p r ~ P l = c 3 

(B3) 


Solutions for these three equations have the following form: 


2 2 4, 


Pr = 



( 2 2) 
^l~ c 2j 

k-p„)+c>„ + c; 

~c\-cl(P' + P u )\ 


[4-4) 


(B4) 


Pr = 


C 2 C 3 + 4 P r ± 2c l c 2 c 3 Jf c l + c 2 ~ P r> “ 4 + 4 P u + 4 [~4 + C V P r + P u\ 


c \ + c 2 


(B5) 


P 


L 


■ c 2 c 3 + 4 P r ±2c \ c 2 c 3^ 

2 1 
c 3 + 1 


-4) 

\( P r- 


■ C 3-C^(P r + P„/ 


r 2 

-4) 

2 


(B 6 ) 


These are potential solutions for equations (Bl), (B2), and (B3), respectively. It is important to realize 
that physically realistic solutions may not exist for certain values of the parameters, c^ 9 c 2 ,c^, and the 
pressures. 
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